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Abstract. We study the difference between on site Hubbard and long range Coulomb repulsions for two 
interacting particles in a disordered chain. While Hubbard repulsion can only yield weak critical chaos 
with intermediate spectral statistics, Coulomb repulsion can drive the two particle system to quantum 
chaos with Wigner-Dyson spectral statistics. For intermediate strengths U of the two repulsions in one 
dimension, there is a crossover regime where derealization and spectral rigidity are maximum, whereas the 
limits of weak and strong U are characterized by a stronger localization and uncorrelated energy levels. 



PACS. 05.45.+b Theory and models of chaotic systems 
systems - 73.20.Jc Derealization processes 



73.10.-w Theories and models of many electron 



1 Introduction 

In low dimensions (d < 2) disorder always yields Q a finite 
localization length L\ when the particles do not interact 
and there is no spin-orbit scattering. When one wants to 
study the role of electron electron interaction, a first issue 
is to know what kind of interaction is appropriate. When 
there are many carriers inside a large length L\ (large 
density and weak disorder), it looks reasonable to assume 
weakly interacting Landau quasi-particles and to take the 
usual short range screened Coulomb repulsion. But, for 
, low carrier densities (10 10 -10 carriers per cm 2 is nowa- 
days achieved in two dimensional heterostructures) the 
screening of the charges is somewhat problematic and one 
may find safer to consider bare long range Coulomb re- 
pulsion. One has in this case a system having charge crys- 
tallization as a natural limit when kinetic energy becomes 
negligible compared to Coulomb energy. The range of the 
interaction can also be varied by metallic gates located 
in the vicinity of the electron gas, as it is often done for 
having a tunable carrier density. This gives us the motiva- 
tion to study the difference between on site Hubbard like 
repulsion and long range Coulomb repulsion in a simple 
limit: two electrons in a disordered chain. Since on-site in- 
teraction plays a role only if the orbital part of the wave 
function is symmetric, we restrict our study to the case 
with opposite spins. 

The problem of two interacting particles (TIP) in one 
dimension has been mainly studied with on site Hubbard 
repulsion of strength U . As proposed by Shepelyansky || 
it has been numerically proven [|,||J|,01 that interaction 
delocalizes a certain number of TIP states over a length 



Li » L\. For a "contact" interaction as Hubbard re- 
pulsion, the TIP system exhibits remarkable properties 
[|||Jl^,|ll|. The mixing of the one body states inside a 
scale L = L\ and the associated derealization effect for 
sizes L > L\ is maximum for U ss U c , where U c is the 
fixed point of a duality transformation || mapping the 
weak U/t limit onto the weak t/U limit, t being the ki- 
netic energy scale. For the one body problem, the spectral 
statistics contains important information: an Anderson in- 
sulator has uncorrelated levels (Poisson statistics) whereas 
a disordered metal displays Wigner-Dyson rigidity charac- 
terizing quantum chaos. At the mobility edge, lies a scale 
invariant critical statistics |T^| , which exhibits a weaker 
spectral rigidity associated to weak critical chaos. For the 
TIP system with Hubbard interaction, the spectrum is 
Poissonian when U —> and U — > oo and becomes [^| more 
rigid when U ~ U c . However, the maximum possible rigid- 
ity does not correspond to Wigner-Dyson rigidity, but to 
an intermediate rigidity analogous to those characterizing 
the one body spectrum at a mobility edge (critical statis- 
tics). It was noticed in Ref. JlO| that those intermediate 
statistics are also related to very slow interaction induced 
TIP diffusion at scales L\ > L > £2- We show in this 
study that, in contrast to Hubbard repulsion, Coulomb re- 
pulsion can drive the TIP one dimensional system to full 
quantum ergodicity with Wigner-Dyson statistics. From a 
statistical study of the interaction matrix elements cou- 
pling two free particle (2FP) states (i.e. the TIP eigen- 
states at U = 0), one finds that Coulomb repulsion mainly 
favors hopping terms between 2FP states nearby in en- 
ergy, with energy separation of the order of the TIP level 
spacing A2 oc L~ 2 . Hubbard repulsion mainly induces || 
hopping terms between 2FP states separated by a larger 
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energy A^^ > A 2 , and the measure of the coupled 2FP 
states is multifractal ||. However, the generic behavior 
of a TIP system with either Hubbard or Coulomb repul- 
sions can be summarized by three regimes: a free particle 
limit dominated by Anderson localization; a large inter- 
action limit dominated again by Anderson localization for 
Hubbard and by charge crystallization for Coulomb;and 
between those two Poissonian limits lies an intermedi- 
ate regime characterized by a maximum mixing of the 
one particle states and a maximum derealization effect. 
The intermediate regime in both cases is located around 
the interaction strengths U c for which the TIP system 
has participation ratios of same order in both preferential 
eigenbases characterizing the weak and strong interaction 
limits. 



2 TIP Hamiltonian 

The TIP Hamiltonian Ti is given by the sum of two terms: 
the first TL, gives the kinetic energy (parameter t) and 
the random potentials (parameter W) in which the two 
particles can move, 

H, = -tJ2 4c j + W^Vim. (1) 

Vi is randomly taken in the interval [—5,5], cf creates 
a particle on the site i and n, = c^q. The second term 
U is the two body repulsion, which can be either on site 
Hubbard repulsion: 

U = V S £ j n l {n l -\) (2) 

i 

or long range Coulomb repulsion: 

w = t/5>K-i ) + | £ (3) 

i i,3=l 1 1 

|i-j|<L/2 

The convention in this work is that two particles at the 
same site cost an en ergy 2XJ (and not U as assumed in 
previous refs. |^,||,^o|JlJ|). An additional cost of energy 
U /p has to be paid by two particles separated by a distance 
p with L/2 > p > 1 when there is Coulomb repulsion. The 
boundary conditions (BCs) are taken periodic. 

Let us denote (e a ,ip a ) and (E a p,ijj a p) the eigenener- 
gies and eigenfunctions of the one particle state |a) and 
of the 2FP state \a) = \a(3) respectively. One has 
E aj3 = e a +€f3. The 2FP level spacing is A 2 « 2B/L(L+1) 
where the band width B « 8t + 2W. The one particle lo- 
calization length L\ is defined from the weak disorder for- 
mula L\ — 100/ VF 2 . Hereafter, the energies will be given 
in units of the kinetic energy hopping term t restricted to 
nearest neighbors. 
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Fig. 1. Density of states p 2 {E) for L = Li = 100. (a) 2FP 
states (U = 0), insert: density around E « (normalized to 1); 
(b) Coulomb ((7 = 5); (c) Hubbard (U = 35) The main band 
of L(L — l)/2 states remains centered around E = together 
with L molecular states at E w 2(7. Coulomb repulsion at (d) 
U = 70 and (e) U = 500. Each subband is centered at E ps U/d 
with d = 0, . . . , L/2. (f) Integrated density showing that each 
subband has L states for Coulomb repulsion (U = 500). 

3 TIP density of states 

When one compares the two repulsions, a first differ- 
ence appears in the density of states 02(E). In the limit 
U — ► 00, Hubbard repulsion splits M the TIP band in 
two parts: a small band of L "molecular states" of high 
energy w U + 2Wvi corresponding to two electrons local- 
ized on the same site i and a main band of L(L — l)/2 
"hard core boson" states which remain at the same small 
energies for U — > 00 and U — > 0. The "hard core bo- 
son" states are given by the resymmetrization |t3| of 
Slater determinants corresponding to electrons in the one 
body state \a > and f3 > respectively. Those states do 
not feel on site interaction, are not coupled to one an- 
other and become decoupled from the molecular states of 
much larger energies when U — > 00. When one takes rigid 
BCs, the resymmetrization is simple and one has exactly 
E a p(U — > 0) = E a p(U — * 00) for a ^ (3. For two electrons 
in a ring enclosing a flux <f> the resymmetrization is more 
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subtle and E a/3 (<j) + <j> /2, U -> 0) = E a p(<p, U -> oo). For 
periodic BCs, E a @ (a ^ (3) goes to the corresponding 2FP 
eigenenergy with anti periodic BCs. In contrast to Hub- 
bard repulsion where the majority of the TIP energy lev- 
els does not feel the interaction when U — > oo, excepted L 
"molecular" states, Coulomb repulsion eventually crystal- 
lizes all the TIP states as two particle "molecules" . When 
L is even, the sizes of the molecules are d = 0, . . . , L/2 
and the spectrum is split in L/2 subbands of L states 
(d L/2) and one subband of L/2 states (d = L/2), each 
of them centered around an energy U/d. When N is odd, 
one has (L+l)/2 subbbands of L states. Without disorder, 
each subband shrinks onto a single L-fold degenerate state 
obtained from successive translations of the "molecules" 
by one lattice spacing, a degeneracy which is broken by the 
random potentials. The different densities pi(E) induced 
by the two repulsions are illustrated in Fig. |l|, for a chain 
of size L = L\ = 100 and various interaction strengths U. 



4 Crossover between two preferential 
eigenbases 

When one turns on Hubbard repulsion, it has been de- 
tailed in Ref. || how the TIP system goes from the 2FP 
basis towards the "hard core boson basis" when U — ► oo. 
The interaction threshold U c = (24) 1 / 4 t/2 was defined 
(for energies near the band center) as the fixed point of 
the duality transformation mapping the distribution of the 
interaction matrix elements oc U/t which couple the 2FP 
states onto the distribution of the kinetic energy matrix 
elements oc t/U which couple the hard core boson states. 
At U m U c , the 2FP basis ceases to be preferential com- 
pared to the hard core boson basis. Being unable to extend 
this duality argument for Coulomb interaction, we study 
the participation ratios PRo of the TIP wavefunctions \^) 
onto the 2FP eigenbasis and PRoo onto the basis built 
out from symmetrized products of site orbitals (site ba- 
sis). This later basis describes the correlated "molecules" 
created when U — > oo. 



Pi?o = (^|(a/3|^)| 4 )- 1 



a/3 



PR 00 = C£\(ij\*)\T 1 



(4) 
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Fig. 2. Participation ratios PR (o) and PR^ (•) of the TIP 

eigenstates of energies ps 2U/L onto the U = and U — oo 
eigenbases respectively. L = L\ = 50. 



coupling the 2FP states. The 2FP wave function have com- 
ponents \P ai p{n, m) on the sites \nm > given by 



< aP\nm >= 
We denote 



ip a (n)ipf3(m) + ii) a {m)ipp{n) 
V2 



?5(0) = E CW^W^tW^W (5) 



n=l 



and 



perm 



QZiite) = E 1 y 

n=l ' y ' 

where perm means the terms obtained after permuting 
(a <-> /?), (7 <-> 6, ) and (a <-> /3, 7 <-> S) for p 7^ 0. 

In the 2FP eigenbasis, the interaction matrix elements 

are 

(a/3\U\76)=WQl S p(0) = UH$ (7) 
for Hubbard repulsion, and 



This allows us to extend for Coulomb repulsion the con- 
cept of a crossover threshold U c where the U = eigen- 
basis ceases to be preferential compared to the U = 00 
eigenbasis. As shown in Fig. |[ one has U c w 120 when 
L = L\ = 50 for Coulomb repulsion. 



5 Interaction matrix elements in the 2FP 
basis 

Before discussing the TIP spectral statistics, it is useful 
to study the structure of the interaction matrix elements 



L/2 

*Y,QZ s (p)) = ucr a 
P =i 



■7 5 



(8) 



for Coulomb repulsion. 

In the absence of disorder and with periodic BCs, the 
one body states are plane waves tpa{n) = (exp ik a n) / \f~L 
and the interaction matrix elements l^jt only couple 2FP 



states of same momentum K = 
For Hubbard, one has when L\ 



"a/3 

K a p = k a 



K. 



H 



7<5 
a/3 



7<5 * 



(9) 
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Fig. 3. Magnitude of the absolute value of the interaction 
matrix elements coupling a given state \aa) (with e a ~ 0) to 
L 2 states |7<5) for a typical sample with L = 2Li = 100. The 
L 2 absolute values are given in the plan (7, 8). The one particle 
states are ordered by increasing values of the energy. The white 
(black) points correspond to maxima (minima) with a linear 
scale of graduation. Top figure: Hubbard H2a\ Bottom figure: 
Coulomb C2 
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The interaction matrix has a block diagonal form. If L is 
odd, one has L blocks of size N s = (L + l)/2. If N is even, 
one has L/2 blocks of size N s = L/2 and L/2 other blocks 
of size N s = L/2 + 1. The N s TIP eigenenergies E n {K) of 
same momentum K are given by the N s solutions of: 



E 

7<5 



1 



E n (K)-E 



7<5 



L 

W 



(10) 



where E 1 $ = 2 cos fc 7 + 2 cos kg. The E n (K) alternate with 
the 2FP energies E 7 s of same momentum. The TIP spec- 
trum is the uncorrelated sum of L such series of different 
momenta K. 

For Coulomb, the previous block diagonal structure is 
preserved when L\ — > 00, but each block becomes more 



4e-03 r- 

3e-03 - 
<|Q(AE)|> 

2e-03 - 
1e-03 - 





10 



100 



1000 



Fig. 4. Disorder average amplitudes < |i?2al > (dashed line) 
and < |C2q| > (continuous line) of the matrix elements cou- 
pling 2FP states separated by an energy AE = \E aa — E 1 s\ 
(in units of A 2 ) for L = 2L\ = 100. 



complex: 



-j-$K al3 ,K., s X 



(11) 



(1 + ^ cxp-^fa - fc 7 )p] + ^ exp -i[(fc/3 - ks)p] 



P =i 



L/2 



L/2 



2p 



+ y-y cxp-z[(fc a - fc 7 )p] + ^ evp-i[(k a - k s )p] ^ 



P =i 



P =i 



2p 



In the presence of a random potential, the absolute val- 
ues of the matrix elements can be given using a linearly 
graduated grey scale in the plane (7, 5) for a given 2FP 
state \a[3 >, the one body states I7 > and \5 > being or- 
dered by increasing energies. When W is sufficiently small, 
TIP momentum remains almost conserved and one can see 
in the plane (7, 5) a white cross made by the two diago- 
nals if a = /3. The first diagonal corresponds to coupling 
to other states \jS > with I7 >« \S >, the second to cou- 
pling to 2FP states \^/5 > close in energy (E a p w E 7 s)- 
When W is larger, Coulomb and Hubbard give rise to 
a different pattern in the (71$) plane: energy-momentum 
conservation remains partially preserved by Coulomb re- 
pulsion (the second diagonal persists) and is lost by Hub- 
bard repulsion, as shown in Fig. || for L w 2Li = 100. 
2FP states \aa > with e a ss are considered in Fig. ||, 
but similar conclusions can be drawn from arbitrary 2FP 
states I a/3 >. 

To explain why C^i continues to mainly couple |a/3 > 
to states 1 7(5 > nearby in energy, we note that when L\ 
is finite, disorder smears the sharp delta function into a 
broader gaussian peak of width a oc L± . (SK a a,K s ) ~~ > 
exp—(K a p — K 7 s) 2 / (2<J 2 ) and (exp— i(k a — kp)p) — > 
exp— (k a — kp) 2 1 (2a 2 ), one gets respectively 



H 



cx exp - 



[K, 



a/3 



K. 



is) 



2a 2 
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and 

Cg«e»p(- ^3g^g)! ) x (12) 

[2 _i n (fcM_i n( fcM) 

a a 
- ln( — - ^ ) - ln( — — ^ )] 

(T (T 

This explains the behavior shown in Fig. ^ where 
the disorder averaged amplitude of the hopping terms 
are given as a function of the energy difference. The L 
hopping terms between the 2FP states |7<5 > nearby in 
energy (AE < A2L/2) are much smaller for Hubbard 
than for Coulomb. For larger energy separation AE, there 
is a \n{\AE / A2 1) decay which is more pronounced for 
Coulomb than for Hubbard. 



6 TIP diffusion and 2FP lifetime 

Coulomb repulsion couples a density « P2 of 2FP states 
nearby in energy. Hubbard repulsion effectively couples a 
smaller density p 2 < P2, as explained in Ref. ||. Let us 
review three consequences of this difference. 

(i) TIP diffusion: In the first studies |,0Jl| of the 
TIP problem, a density p2{Lx) of 2FP states coupled by 
the interaction was assumed for L « L\. Under this as- 
sumption, it was predicted that the TIP dynamics should 
exhibit interaction assisted diffusion on scales L\ < L < 
L2, the time evolution of the TIP center of mass R2 being 
given fjJI by: 

R 2 (t) « y/D 2 (t)t (13) 

with D2(t) is roughly constant, up to log(t) corrections 
In Ref. a much slower propagation R2 oc logi 
was observed for Hubbard repulsion, attributed to the 
weak density p 2 (Lx) of effectively coupled states. One 
expects that the original prediction will be at least par- 
tially restored for Coulomb repulsion. 

(ii) TIP localization: The interaction assisted propaga- 
tion stops at a scale L2 characteristic of TIP localization. 
Assuming | < af3\U\jS > | 2 oc U 2 /L\ for L w Lx, the 
enhancement factor L2/ L\ was originally given by the es- 
timate: 

^(x\<a(3\U\ 1 S>\ 2 p 2 (L 1 )ocL 1 (14) 
Lx 

It was pointed out in Ref. S that, since p 2 (Lx) oc L\ — > 
p2 (Li) oc L{- 75 in the presence of Hubbard repulsion, the 
enhancement factor should be weaker (L 2 /Lx oc \fL~x), a 
prediction confirmed by numerical calculations. It is likely 
that L2/L1 oc Li will be a better estimate for Coulomb 
repulsion. 

(iii) 2FP lifetime: The inverse lifetime (r af3 ) of a 2FP 
state I a/3) is given by Fermi Golden rule: 

r a0 oc I < a0\U\j8 > \ 2 6(E af} - E lS ) (15) 

7<5 
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L, 

Fig. 5. Inverse lifetime r aa (Li) of a 2FP state \aa > with 
e a ~ for Hubbard (o) and Coulomb (o) repulsion respectively. 
L = 420. The continuous (dashed) line corresponds to r aa cx 
L- 1 * 5 (ccL- 1 - 05 ). 

If the density of coupled states is p 2 (Li) « L\ , a reason- 
able estimate for Coulomb repulsion, one gets 

r aP (L = L x ) oc I < a/3\U\j6 > | 2 p 2 (ii) cx L^ 1 . 

In contrast, the lifetime should be longer for Hubbard re- 
pulsion: 

r aP (L = Lx) oc I < af3\U\l$ > \ 2 p e 2 "(Li) oc L^ 1 - 65 

since the density j|) of coupled states (by the second mo- 
ment of I < a/3\U\-yS > |) is p e 2 ff oc i} 35 for a 2FP state 
\aa >. 

We have checked this prediction. In Fig. ||, the inverse 
lifetime calculated in a chain of length L — 420 is shown 
as a function of L\. When Lx < L, one has r aa oc i^ 1 ' 65 
for Hubbard repulsion while r aa oc i^ 1 ' 05 for Coulomb 
repulsion, close to a simple decay oc L± . When L becomes 
larger than Lx, the lifetime becomes L-independent for 
Hubbard, and may continue to weakly decay as a function 
of L for Coulomb. Fig. || shows us that this possible decay 
remains negligible. 



7 Two particle spectral statistics 

We now study how the spectral statistics depend on the 
interaction strength U for the two repulsions, using the 
distribution P(s) of energy spacings between consecutive 
levels and the variance E2 (E) of the number of levels in- 
side an energy window of width E. We consider energy 
levels in the bulk of the low energy sub-band: E ss for 
Hubbard repulsion and E ss 2U/L for Coulomb repulsion 
(see Fig. |l|). We take L = Lx for having the largest inter- 
action matrix elements, and hence the maximum mixing 
of the 2FP states. 
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Fig. 6. Spacing distribution P(s) for Hubbard (•) and 
Coulomb (o) repulsions at U = 70 and L — L\ = 50. The 
symbols (o) correspond to a medium range interaction trun- 
cated to p = 1. The levels are taken around E as 2U/L for 
Coulomb repulsion and around E as otherwise. The continu- 
ous lines correspond to Pp{s) and Pw{s) respectively. Insert: 
Parameter r\ around E as as a function of the range p of the 
interaction. 

After unfolding the spectra, one expects 

P(a)«iV(*)« yexp(~5 a ) (16) 

and 

S 2 (E) as Sf{E) as -1(^(2^) + 7 + i - (17) 

7T^ 8 

(7 being the Eulcr constant) for correlated levels having 
Wigner-Dyson statistics, whereas one should have Poisson 
statistics with 

P P {s) = exp(-s) (18) 

and 

IJP(E) = E (19) 

for uncorrelated levels. 

The TIP spectra are well described by Wigner-Dyson 
statistics for an intermediate Coulomb repulsion. In Fig. || 
and in Fig. 0, one can see indeed that the level spacing dis- 
tribution and the number variance S 2 are well described 
by the Wigner surmise Pw(s) and (E) respectively 
near the intermediate interaction strength U c as 120 for 
which the TIP system does not have a preferential eigen- 
basis (see Fig. ||. As shown by S 2 (E) (Fig. 0), Wigner- 
Dyson spectral rigidity is established over an energy inter- 
val containing a few levels. In contrast, Hubbard repulsion 
can only yield || at the corresponding U c (as 1) 

P SP as 4sexp-(2s) (20) 

and 

S 2 ~ 0.16 + 0.41P (21) 




AE/A 2 



Fig. 7. Number variance £2 as a function of the energy in- 
terval E in units of TIP level spacing A2 (same symbols as in 
Fig. P excepted the filled diamonds (Coulomb at U = 100 and 
L = Li = 100). 

respectively. 

To study how the spacing distribution depends on U, 
we use the spectral parameter 77 defined by 

x J^ds[P{s) -Pw(s)] , x 

^ds[P P ( S )~P w ( S )} 

with b = 0.4729. For {7 = 0, the consecutive levels are 
essentially uncorrelated and rj at 1. The curves rj{U) given 
in Fig. H show us a striking difference between the spectral 
statistics yielded by the two repulsions. 

For L = L\ = 100, let us estimate the interaction 
threshold U* for which the interaction matrix elements 
coupling consecutive 2FP states becomes of the order of 
their energy separation A 2 . p 2 ~ 0.25 (density normalized 
to 1) gives A 2 « (p 2 L(L + l)/2)~ 1 « 4/5050 (see insert 
of Fig. for the energy spacing between consecutive 2FP 
levels around E = 0. The off-diagonal interaction matrix 
elements Q = H^p (Hubbard) or Q = C^t (Coulomb) 
coupling consecutive 2FP levels are normally distributed, 
as shown in Fig. |^. The root mean square of Q is small 
for Hubbard (w 0.0038) and larger for Coulomb (as 0.01). 
This gives U* c as 0.08 for Coulomb and U* H as 0.2 for 
Hubbard. 

When U increases but remains lower than U*, there 
is first a perturbative regime discussed in Refs. [p^ , pT| 
where the interaction yields Rabi oscillations between con- 
secutive 2FP states at a frequency given by the absolute 
value of the coupling interaction matrix element. More- 
over, the energy range Ejj under which one has level re- 
pulsion is given [^6j by this Rabi frequency. When U be- 
comes equal to U*, one has a transition from this per- 
turbative regime towards an effective Fermi golden rule 
decay of the 2FP states and the characteristic ran ge En 
over which Wigner-Dyson rigidity occurs becomes [16ln7l 
proportional to the square of the amplitude of the coupling 
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Fig. 8. Distribution P(Q) of the interaction matrix elements 
Q — H^fg (Hubbard: circle) or Q — (Coulomb: triangle) 
coupling nearest neighbor 2FP states (energy separation ~ A?) 
for L = L\. The distributions are fitted with gaussian curves of 
variance a 1 = 0.01 2 (Coulomb, dashed line) and a 2 = 0.0038 2 
(Hubbard, continuous line). 

1.00 | — — — — ^ — 

<> / 




Fig. 9. Spectral parameter r\ as a function of U for L — L\. 
Levels around E ~ 2U/L for Coulomb repulsion with L = 
100(»), 50(o) and 65(A). Levels around E « for Hubbard 
repulsion (+) with L = 100. 

matrix elements: Ejj oc U when U < U* and Ejj cx U 2 
when U > U* . U = U* is the interaction threshold where 
the spectral statistics are intermediate between Poisson 
and Wiener (77 w 0.39 when P(s) w Psp(s)). Looking 
at Fig. || one can see that 77 decreases as a function of 
U down to the characteristic value 77* sa 0.39 reached 
when U ~ C/^ w 0.08 for Coulomb repulsion, and when 
U ~ Uff rs 0.2 for Hubbard repulsion. 

For Coulomb repulsion, 77 continues to decrease when 
U > Uq down to 77 = with the slower U dependence 
characteristic of the effective golden rule decay of the 2FP 



states. The Wigner-Dyson distribution P\y(s) is fully es- 
tablished at U rs U c , i.e. when the TIP system is exactly 
as far from the U = eigenbasis than from the U = 00 
eigenbasis. When U > U c , integrability is slowly restored 
and 77 — > 1 as U — > 00. 

For Hubbard repulsion, the spectral rigidity does not 
continue to increase above U^, but saturates to the inter- 
mediate critical rigidity characterized by Eq. ^o| and Eq. 
[H]for P(s) and ^(-E) respectively. Above the fixed point 
U c ~ 1 of the duality transformation, the TIP system 
becomes closer to the U = 00 eigenbasis and the levels 
become statistically uncorrelated. This critical Hubbard 
regime is a complicated issue where the multifractal char- 
acter of the interaction matrix described in ref. [g] is very 
likely relevant. However, one can do the following remark. 
When L\ — * 00, the interaction matrix is block diagonal, 
a block corresponding to a pair momentum K. For Hub- 
bard, the N s TIP level E n (K) of momentum K are located 
near the 2FP levels of same momentum when 4XJ /L — > or 
AU/L — > 00. One can assume that they should be near the 
middle of consecutive 2FP levels of same K for U ~ U c . 
The distribution P(s) of levels located in the middle of 
levels with spacing distribution Pp{s) is the semi- Poisson 
distribution Psp(s). If the sequence of 2FP levels of mo- 
mentum K were randomly distributed, the spacing distri- 
bution of the N s TIP levels should be given by Psp(s) 
without disorder. On the contrary, the interaction matrix 
being more random for Coulomb, one can expect that the 
L series of N s levels of momentum K will be driven to- 
wards Wigner-Dyson statistics as U increases. However, 
though the argument may give hints for the existence of 
the Hubbard Pgp(S), it does not explains the behavior of 
S%{E). The breakdown of momentum conservation by the 
disorder, and the associated mixing of the L independent 
series of N s levels characterizing the clean limit plays a 
complex role. 

The insert of Fig. || shows how the spectrum becomes 
more rigid at U — 70 when the range p of the interaction 
is increased. S^iE) displays a similar information in Fig. 

& 



8 Quantum melting for intermediate 
Coulomb repulsions 

The intermediate Wigner-Dyson regime yielded by Coul- 
omb repulsion corresponds, inside a scale L%, to a complete 
melting of the localized 2FP states previous to crystalliza- 
tion. To show this, we introduce two parameters 7 and 
£. For a TIP wavefunction |^), we calculate the density 

Pi = (t/?|cfc,;|!/?) at site i and the density density corre- 
lation function C(r) = (l/2)^2i Pipi-r- The participation 
ratio £ (i. e. the number of sites occupied by a TIP state) is 
given by 2C(0)~ 1 . The crystallization parameter 7 is given 
by the difference MaxC(r) — MinC(r), where all trans- 
lations r (including r = 0) are considered. If the electron 
density is homogeneous (as for an extended liquid state) 
7 w whereas 7 sa 1 if the two charges are mainly located 
on two different sites, a situation occuring when U — > 00 
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p8| of one dimensional spinless fermions at half filling. In 
two dimensions, there are experimental evidences that 
localization may disappear for intermediate Coulomb en- 
0.8 ergy to Fermi energy ratios r s . A new quantum regime 
has been observed |19| ] for the ground state of two dimen- 
A jj ^/L sional spinless fermions at intermediate factors r s where 

j the metallic phase is observed. The nature of the states, 
J ~ °' 6 when the system is far from the free particle (Fermi glass) 
| £ or the Coulomb (Wigner crystal) bases, and the associ- 
, ated transport mechanism, remain to be understood. This 
a • A simple one dimensional study draws our attention to the 
A A ' important role played by the range of electron-electron 
g a i repulsions. To know how is or is not screened Coulomb 
^ repulsion in low density electron systems is then an im- 
i „ „ portant issue. 
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Fig. 10. Crystallization parameter 7 for L = L\ = 50(A) 
and 100(o) and the corresponding numbers of occupied sites £ 
(filled symbols) for Coulomb repulsion for TIP levels around 
E = 2U/L. 

(Coulomb "molecule") and when C/ — > and L\ — > (two 
electrons located in two minima of the potential). The 
variations of £ and 7 are given in Fig. [To] when Coulomb 
repulsion U increases. The curves y(J7) and £(U) are cor- 
related to the curve r](U) shown in Fig. ||. The curve £(U) 
shows us that the TIP states occupy a fraction of the chain 
without repulsion before being uniformly spread over a 
scale L\ at interaction strength U ~ U c for which there is 
a crossing of the two curves PR(U) given in Fig. |^. For 
L = L\, one has an interaction induced quantum melting 
of the non interacting glass, yielding quantum ergodic- 
ity with " more liquid" and extended wavefunctions and 
Wigner-Dyson spectral statistics. 



9 Conclusions 

We have studied one of the simplest problems where quan- 
tum localization and two body interaction are in com- 
petition. We have seen that the range of the interaction 
makes important differences. One of them is revealed by 
the study of the spectral statistics, providing an intriguing 
puzzle for quantum ergodicity: the Wigner-Dyson statis- 
tics shrinking to intermediate statistics when the two body 
repulsion becomes local. However, the generic behavior of 
the TIP system can be summarized by three regimes, inde- 
pendently of the interaction range. There is the free parti- 
cle limit dominated by quantum localization, the Coulomb 
limit dominated by the pinning of a correlated system 
of charges ("Coulomb molecule, Wigner crystal") or by 
quantum localization again (Hubbard repulsion) . Between 
these two limits, there is an intermediate regime where one 
has a maximum mixing of the one body states, making the 
states more extended and the spectrum more rigid. In one 
dimension, this yields a partial derealization effect [hi 
could be large but remains finite). Similar conclusions have 
been reached from a study of the many body ground state 
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